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Abstract. The Next-To-Minimal-Supersymmetric extension of the Standard Model (NMSSM) has 
been in the focus of extensive studies in the past two decades. In anticipation of the LHC era, 
O^l the interest in automatized tools that can calculate collider signatures has grown. We present the 

implementation of the NMSSM into the event generator WHIZARD. In addition to a brief review 
of the implementation, we discuss the testing and validation procedure. Phenomenological studies 
\ will not be presented here. 
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'clj THE MULTI-PURPOSE EVENT GENERATOR WHIZARD 

WHIZARD [1] is a multi-purpose Monte-Carlo event generator for the Standard Model 
(SM) and beyond (BSM). It uses a multi-channel adaptive phase space integration based 
on the VAMP algorithm [|20, which is called by a highly efficient phase space grid de- 
composition relying on heuristics of the underlying resonance structure of field theoretic 
scattering amplitudes. The matrix element are delivered by the generator O'Mega [H] 
which avoids in an optimal way all redundancies of tree-level scattering amplitudes and 
^ ! combines gauge-invariant substructures. WHIZARD was the first generator for full ma- 

trix elements for the MSSM [4] and has been validated by many groups and projects both 
for supersymetric models and alternative BSM models [5]. Among the available mod- 
els are half a dozen variants of Little Higgs models, extra-dimensional models, several 
(extended) supersymmetric models, general anomalous gauge, top and Higgs couplings, 
general unitarization models for electroweak scattering processes, noncommutative ver- 
^ . sions of the Standard Model etc. 



WHIZARD can be accessed via the Hepforge web site: |http : / /projects . hepf orge . org/whize 
The alpha version of release 2.0.0 has been finished just before the SUSY conference, 
including the first implementation of the Next-to-Minimal Standard Model (NMSSM) in 
a multi-particle event generator, which is the main topic here. We briefly want to review 
the new features of version 2.0.0, which brings a major improvement over the older 
version 1 . WHIZARD has been completely restructured in a fully object-oriented way. 
It allows for an arbitrary setup of structure functions, comprises a fully flexible script 
language for defining scales, cuts, analyses variables as arbitrary functions of kinematic 
variables and particle lists, it contains a rudimentary own parton shower (which will 
be further developed during the development of WHIZARD 2 and connected with 
a module for simulating multiple interactions which is right now being developed). 
WHIZARD is able to write out standard event formats like LHA, LHEF, STDHEP, 
HepMC etc. It comes with its own graphics analysis tool, while an interface to ROOT 



will also be available in WHIZARD 2. Also, an interface to DELPHES [6] for a fast 
detector simulation is under way. 

Although, WHIZARD is able to deliver matrix elements for 2 — > 16 processes (or 
higher multiplicities) and integrate processes with at least up to ten particles in the final 
states, this might not be sufficient for BSM models with a discrete parity symmetry as 
suggested by the existence of dark matter. Hence, WHIZARD 2 is also able to integrate 
and simulate long cascade decay chains with full spin correlations. A further improve- 
ment is a full-fledged interface to FeynRules [0] for a more automatized inclusion of 
new physics models for both WHIZARD versions 1 and 2. 



THE NMSSM 

The Next-to-Minimal Supersymmetric Standard Model (NMSSM) is a viable extension 
of the MSSM in the Higgs sector of the model (see e.g. the references in [8]). The 
main motivation for the NMSSM is a possible solution to the so-called /i problem, 
the fact, that the supersymmetric parameter /i in the superpotential is dimensionful and 
should in principle be of the order of the SUSY breaking scale. A working electroweak 
symmetry breaking demands this parameter, however, in the ball park of a few hundred 
GeV. The NMSSM tackles this problem by enlarging the particle spectrum by a single 
left-chiral superfield 5, being a singlet under the SM gauge group. Their are two general 
options to generate a quartic potential term for the scalar part of S: either a D term of 
an additional U(l) guage symmetry (like e.g. in the PSSSM [9]) or by a F term from a 
cubic superpotential term (conventional NMSSM). For the conventions, we follow those 
of the SUSY Les Houches Accord 2 [10] (for a recent review on NMSSM conventions 
cf. also [8|]). This is a quite general approach, but neglecting possible CP, 7?-parity, or 
flavor violation. The suporpotential of the NMSSM is given by 

W N mssm = Wmssm - e ab XSH^ + ^S 3 + ^i'S 2 + , (1) 

where we neglect an explicit jl' term and the Fayet-Iliopoulos term for the singlet 
superfield. The most general soft-breaking terms for the NMSSM 

V S0 ft = V 2 ,MSSM + V3MSSM + ml\S\ 2 + (-e ab ?iA x SHtH^ + ^KA K S 3 + m'iS 2 + ^sS+h.c.) , 

(2) 

using the conventions from [ | 1 QT | . 

The field content of the NMSSM is almost the same as for the MSSM, except for an 
additional scalar and pseudoscalar Higgs boson, denoted by H® and A®, as well as a fifth 
neutralino, x<j coming from the additional singlino component. 

IMPLEMENTATION AND VALIDATION 

As mentioned in the section about the model before, we are sticking to the SLHA 2 
conventions. This fixes the uncertainties in signs and phases. As for the MSSM, the 
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FIGURE 1. Comparison between WHIZARD and the FeynRules implementation within Madgraph. 
The green color shows agreement, at least up to the Monte Carlo integration error. We used y/s = 3 TeV 
to be above all thresholds. 



WHIZARD implementation uses explicitly positive masses for charginos and neutrali- 
nos, and puts signs and complex phases instead in the corresponding mixing matrices. 
In contrast to the MSSM, the NMSSM implementation is more flexible than the MSSM 
one (although now thanks to Bjorn Herrmann there is a completely general MSSM im- 
plementation in WHIZARD), as it allows not only for a full CKM matrix but also for a 
left-/right mixing for all generations, but no inter-generational mixing. The generaliza- 
tion to the CP-non-conserving case is easily possible. 

The NMSSM is a whole is an incredibly complicated model with order 6,700 cou- 
plings (including quartic and Goldstone couplings, most of each are fortunately not of 
phenomenological importance). It is mandatory to check an implementation as far as 
possible. For this task, we followed the strategy given in Ollh : unitarity checks for 2 — > 2 
and 2 — > 3 scattering amplitudes have been performed, and we tested Ward- and Slavnov- 
Taylor identities for both gauge and supersymmetry. Another check was to reproduce the 
correct MSSM limit for the NMSSM, namely letting the trilinear singlet Higgs coupling 
and the cubic singlet coupling approach zero, sending the singlet vev to infinity while 
keeping the combination of the latter two fixed to /1: A — > 0, K — > 0, (S) — > 00, (S) A — > /i. 

What is by far one of the most stringent tests is the comparison with an independent 
implementation. For that purpose we are using a FeynRules generated NMSSM model 
file for MadEvent/Madgraph [12] as well as CalcHEP jl3\ \. Tabled] shows the compar- 
ison for t + t~ initial state processes as an example. The corresponding SLHA 2 input 
file is left out here due to reasons of limited space. Our comparison contains far more 
than 600 processes, showing full agreement between the two implementations. 



CONCLUSIONS 



We presented the complete implementation of the NMSSM implementation into the 
multi-purpose event generator WHIZARD. A vast sample of consistency checks and 
comparisons with independent implementations show the correctness of the code. 
WHIZARD is available as a modern high-level tool for NMSSM signal and background 
simulations for the LHC era. A detailed version of the comparison and validation of the 
NMSSM implementations will appear soon 
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